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Abstract 

In this paper, we are concerned with the global existence and optimal rates 
^ ■ of strong solutions for three-dimensional compressible viscoelastic flows. We 

OO . prove the global existence of the strong solutions by the standard energy method 

under the condition that the initial data are close to the constant equilibrium 
ly^ ■ state in i^'^-framework. If additionally the initial data belong to L^, the optimal 

Q>^ ■ convergence rates of the solutions in L^-norm with 2 < p < 6 and optimal 

I convergence rates of their spatial derivatives in L^-norm are obtained. 
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1 Introduction 

In this paper, we are interested in three-dimensional compressible viscoelastic flows 
[3,5,12,23]: 

Pt + div{pu) = 0, (1.1a) 

{pu)t + div{pu ®u) - pAu - (A + p)Vdivu + VP(p) = adiv{pFF'^), (1.16) 

Ft + u-VF = VuF, (1.1c) 

for {t,x) G [0,+oo) X M^. Here p, u eM.^, F e M^^^ (the set of 3 x 3 matrices with 
positive determinants) denote the density, the velocity, and the deformation gradient, 
respectively. The Lame coefficients p and A are satisfied the physical condition: 

p>0, 2/i + 3A>0, 
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which ensures that the operator — /xA — (A + /x)V(iiv is a strongly eUiptic operator. The 
pressure term P{p) is an increasing and convex function of p for p > 0. The symbol 
denotes the Kronecker tensor product, means the transpose matrix of F, and 
the notation u ■ VF is understood to be {u ■ V)F. For system (1.1), the corresponding 
elastic energy is chosen to be the special form of the Hookean hnear elasticity: 

a \ C 

W{F)^-\F\^ + - P{s)ds, a>0, 
2 P Jo 

which, however, does not rcdTicc the essential difficulties for analysis. Indeed, all the 
results we describe here can be generalized to a more general cases. 

In this paper, we investigate the Cauchy problem of system (1.1) with the initial 
condition: 

{p,u,F)\t=o^{po{x),uo{x),Fo{x)), xeR\ (1.2) 

We also assume that 

div{pF^) = 0, F^''{Q)ViF'^{Q) = F^^{Q)ViF'\Q). (1.3) 

It is standard that the condition (1.3) is preserved by the flow, which has been proved 
in [7,21]. 

For the incompressible viscoelastic flows and related models, there are many im- 
portant progress on classical solutions, refer to [1,2,9,13,16] and references therein. On 
the other hand, the global existence of weak solutions to the incompressible viscoelas- 
tic flows with large initial data is still an outstanding open question, although there 
are some progress in that direction [15,17,18]. For the compressible viscoelastic flows, 
to our knowledge, there are few results on the dynamics of global solutions to com- 
pressible viscoelastic flows, especially on the large time behavior. The local existence 
of multi-dimensional strong solution was obtained in [6], and the global existence of 
strong solution with the lowest regularity was shown in [7,21]. For the initial boundary 
value problem, global in time solution was proved to exist uniquely near the equilibrium 
state in [8,22]. 

In this paper, we firstly study the optimal time-decay rate of the global strong 
solutions to the Cauchy problem (1.1)-(1.2). To be more precise, the main purpose of 
this paper is to study the existence and uniqueness of global strong solutions and in 
particular the asymptotic behavior on the Cauchy problem of compressible viscoelastic 
fiows. We prove the global existence of strong solutions by the standard energy method 
in spirit of Matsumura and Nishida [19,20]. In order to obtain the linear time-decay 
estimates, we need to analysis the properties of the semigroup, as in [10,11,14,24]. 
Unfortunately, it seems untractable, since the system (1.1) has thirteen equations. To 
overcome this difficulty, we take Hodge decomposition of the linear system, then it 
becomes two similar systems, each of those only involves two variables, which makes 
us be able to obtain the optimal time-decay estimates. 

Our main results are formulated in the following theorem: 
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Theorem 1.1. Assume that the initial value {po — 1,Uo,Fq — I) e H'^{M.^) satisfies the 
constraints (1.3), then there exists a constant Sq such that if 

\ipo-l,Uo,Fo-I)\H2<So, (1.4) 

then there exists a unique globally strong solution (p, u, F) of the Cauchy problem 
(1.1) - (1.2) such that for any t G [0, oo), 

\{p-l,u,F-I){-,t)\j,,+ [ \V{p,F)\l, + \Vu\l,<C\{po-l,Uo,Fo-I)\m. (1.5) 

Moreover, if (po - 1, uq, Fq - I) e L^{R^), then 

\(p-l,u,F-I)(t)\Lr<C(l + t)-"^^'-'^\ V pe[2,6], (1.6) 
|V(p -l,u,F- /)(t)|^i < C{1 + t)-t. (1.7) 

Finally, denote 

{go, mo, J^o) = {po - 1, po^o, Po-fo - I) 
and assume that the Fourier transform {go,rno,J^o) satisfies 

|^'~o|>co, |mo|<|er, \J'o-^o\<\C\\ /or < 1^1 «1, (1.8) 

where Co and rj are two positive constants. Then we also have the lower bound time 
decay rate as 

|(p-l)(t)U. >ci(l + t)-i, (1.9) 
|«(t)|L2 >ci(l+t)-t, (1.10) 
|(F-/)(t)|^. >ci(l + t)-t, (1.11) 
where ci is a positive constant independent of time. 

Notaions. We denote by L^, VF'"'^ the usual Lebesgue and Sobolev spaces on 
and — VF"*'^, with norms \ • \lp, | • Iw^^'P and \ ■ \h'^ respectively. For the sake of 
conciseness, we do not distinguish functional space when scalar- valued or vector-valued 
functions are involved. We denote V = = (81,82,83), where 8i = 8^^, = 8i and 
put 8lf = V'/ = V(V'~^/). We assume C be a positive generic constant throughout 
this paper that may vary at different places and the integration domain will be al- 
ways omitted without any ambiguity. Finally, (•, •) denotes the inner- product in L^(R^). 

The rest of this paper is devoted to prove Theorem 1.1. In Section 2, we first refor- 
mulate the system and do some careful a priori estimates for the strong solutions. Then 
the global existence of the strong solutions is established by the standard continuity 
argument. In Section 3 we will derive the decay- in-time estimates for the linearized 
system and use the energy method to derive a Lyapunov-type energy inequality of all 
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the derivatives controlled by the first order derivatives, then we utilize the decay-in- 
time estimates for the linearized system to control the first order derivatives by the 
higher order derivatives. Hence, the optimal decay rates of the global strong solutions 
follow from these two kinds of estimates. In section 4, we establish the lower bound 
time decay rate for the global solution. 

2 Global existence 

2.1. Reformulation 

In this subsection, we first reformulate the system (1.1). Without loss of generality, 
we assume -P'(l) > 0, and denote Xo = (-P'(l))~^ For P > 0, system (1.1) can be 



rewritten as 

Pt + pdivu + uV p — 0, (2.1a) 

1 P'(o) 

ul + u-Vu'- -{pAu' + (A + p)Vidivu) + -^^iP = aF^^V^F'^, (2.16) 

Ft + u-VF^ VuF, (2.1c) 



where we used the condition div{pF'^) — for all t > which ensures that the i-th 
component of the vector div{pFF'^) is 

VjipF'^Fj'') = pF^'^VjF'^ + FikVjipF-'''') = pF^^VjF'^. 

Denote 

n{t, x) = p(xoi, Xox) - 1, v{t, x) = Xo^^(Xo^, Xo^), E(t, x) = F(xlt, Xox) - I, 
then 

rit + div V — f — V ■ Vn, (2.2a) 
vi - pAv' - (A + p)Widiv V + Win - aVjE'^ = g, (2.26) 
Ef-Vv = h-v- VE, (2.2c) 
where ^ 

/ = -nV -v, h = VvE, a = p7^> 

= aE^'VjE^' - -^{pAv' + (A + p)Vdtvv) -vW- ^Jl, - l] V,n. 

l + n \[1 + n)F'[l) J 

Again, without loss of generality, we will assume that a = 1 for the rest of this paper. 

2.2. A priori estimate 

As a classical argument, the global existence of solutions will be obtained by combin- 
ing the local existence result with a priori estimates. Since the local strong solutions can 
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be proven by standard argument of Lax-Milgram theorem and the Schauder-Tychonoff 
fixed-point as [6] whose details we omit, global solutions will follow in a standard 
continuity argument after we establish (1.5) a priori. Therefore, we assume a priori 
that 

\{p-1,u,F-I)\h2<5o<:1, (2.3) 

which is equivalent to 

\{n,v,E)\H2<5<.l. (2.4) 

Here So ^ 5 is small enough. This, together with Soboles's inequality, implies in 
particular that 

\{n,v,E)\Lo. <CS. (2.5) 

This should be kept in mind in the rest of this paper. 

For later use we first estimate the norm of f,g,h. By (2.4), (2.5), together with 
Sobolev's inequality, Hoder's inequality and Moser-type's inequality, we easily deduce 
that 

\{f,h)\L2 < \{n,E)\Loo\Vv\L2 < C5\Vv\l2, 

\V{f,h)\ < \{n,E)\L^\VML'^ + \V{n,E)\L-s\Vv\L. < C6\V\\l2, 

|V2(/,/i)| < C{\{n,E)\Lo.\Vh\L2 + \V\n,E)\L2\Vv\L^) < C6\V\\hi, 

\9\l^ < C(|^;|l3|Vv|l6 + \n\^\V^v\L2 + \{n, E)\Loo\V{n, E)\l2) 
<C5{\VML^ + \V{n,E)\L2) 

|Vc/|l2 < C{\Vv\l--\V^v\l2 + \Vv\ l°°\'Vv\l2 
+ \n\Lo^\V^v\L2 + |Vn|i6|V^f |l3 
+ \{n,E)\Lo.\V\n,E)\i^2 + \V{n,E)\i^e\V{n,E)\Ls) 
< C5{\V^v\H^ + \V\n,E)\L2), 

where we used the fact 

In what follows, a series of lemmas on the energy estimates is given. Firstly the 
energy estimate of lower order for {n, u, E) is obtained in the following lemma. 

Lemma 2.1. Under the priori assumption (2.4), we have 

\j^\{n,v,E)\l,-rC\Vv\' <C8\V{n,E)\l,. (2.6) 

Proof. Multiply (2.2a), (2.2b), (2.2c) by n, v, E respectively and then integrating them 
over R^, we have 

ill(^>^> E)\h + /^|v^^li2 + (// + A)|V • v\l, 

= if -vVn,n) + {g,v) + {h-vVE,E). ^ ' 

The three terms on the right hand side of the above equation can be estimated as 
follows. 
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First, it holds that 

{f — V ■ Vn, n) = (— nV ■ v — v ■ Vn, n) 
— {v ■ Vn, n). 

It follows from Sobolev's inequality. Holder's inequality and (2.4) that 

\{f - V ■Vn,n)\ < |n|i3|^;|i6|Vn|i2 < C|n|i7i|Vu|L2|Vn|i2 

<C5(|Vn|i. + |V^li.)- ^ ^ ^ 

Similar to the proof of (2.8), we have 

\{h-vWE,E)\ < C5{\VE\l2 + \Vv\l2). (2.9) 
For the second term, we have 



\{g\v^)\ = C(|(i?^■^V,i?^^^;^)| + \{^Av\v^)\ + \{^ydtv,v,v 

+ \{v ■ Vv\v^)\ + |(((^:§{^ - l)V,n,^^)|). ^"-^'^ 



As the proof of (2.8), it follows from Sobolev's inequality. Holder's inequahty and (2.4) 
that 

\{E^'\/,E^\v^)\<CS{\WE\l, + \Vv\l.), (2.11) 

\{V ■ Vv,v)\ < \v\Ls\v\Le\Vv\L^ < C\v\Le\Vv\'^ < C5\Vv\l2, (2.12) 

p'(n + ^) 

'^^ (l + n)ni) " - ^I^U^I^I^^I^^U^ ^ C6i\Vn\l, + (2.13) 

|(^A^;\^OI = |(V(^)V^^^^)| + (^V^;^V^;^) 

< C(|Vn|j^i|VT;|L2 + \n\Loo\Vv\L2) (2.14) 

< C5|V^;|i2, 

and similarly, 

71 

\{,^-^\/divv\v') \ < CS\\/v\l2. (2.15) 
Substituting (2.11)-(2.15) into (2.10) gives that the second term is bounded by 

\{g,v)\<CS{\V{n,E)\l, + \Vv\l,). (2.16) 

Hence combining (2.7), (2.8), and (2.16) yields (2.6) since 5 > is sufficiently small. 
This completes the proof of the lemma. 

In the following lemma we give the energy estimate of the higher order for (n, v, E). 

Lemma 2.2. Under the assumption (2.4), we have 

~\V{n,v,E%, + C\V'v\jj^<C5\V{n,E%,. (2.17) 



Proof. Applying V to (2.2a), (2.2b), (2.2c) and multiplying by Vn,Vv,VE respec- 
tively, integrating over M^, we have 

i||V(n, V, E)\l.. + fAV^vll, + (/i + A)|V(V • v)\l, 
= {^{f-v-Wn),Wn) + {Wg,Wv) + {W{h-vWE),WE). ^ ' 

Now let us estimate the right-hand side term by term. First of all, by Holder's inequality 
and Sobolev's inequality, we have 

|(V/,Vn)| + |(V/i,V£;)| < \V{f,g)\L2\V{n,E)\L^ 

< C5\V^v\L2\Vn\L2 

< CS{\V'v\l, + \Vin,E)\l,). 

Next, integrating by parts, we get 

|(V^, V^)| = |(^, V^l < C\g\L2\V'v\L2 < CS{\V'v\l, + \V{n,E)\l,). 

Finally by symmetry, we have 

I {V{v ■ Vn) , Vn) I + I {V{v ■ VE) , VE) \ 
= \{Vv- Vn, Vn) + {v ■ V Vn, Vn) | + | (Vu • \/E, VE) + {v ■ WE, VE) \ 
= I (Vf ■ Vn, Vn) + l{div v, | Vnp) \ + \{Vv ■ VE, VE) + \{div v, \WE\^) \ 
< C\Vv\Lo.\V{n,E)\l, < C6{\V^v\l, + \V{n,E)\l,). 

Substituting these results into (2.18), we conclude 

~\V{n,v,E)\l,+C\V'v\l. < C6{Mn,E)\l, + (2.19) 

Similarly, applying V^ to (2.2a), (2.2b), (2.2c) and multiplying by V^n, V^v, V^-E 
respectively, integrating over M^, wc have 

li\V\n,v, E)\l, + + (^ + A)|V(V2 ■ v)\l, 

= {V^{f-vVn),V^n) + {V^g,V^v) + {V\h-vVE),V^E). ^^■^^> 

To estimate the right-hand side of the above equation, we note, by Holder's inequality 
and Sobolev's inequality, that 

\{V'f,V'n)\ + \{V%V'E)\ < \V\f,h)\r.2\V\n,E)\L. 

<C6{\V'v\l, + \V\n,E)\l,). 

Integrating by parts, we have 

KV^^, VMI = \{^9,^'v)\ < \^9\l2\^'v\l2 < C6{\WMm + \^\n,E)\l,). 
Finally, by symmetry, we have 

\{V\vVn),V^n)\ 
= I ( V^f ■ Vn, V^n) + ( V^; ■ V Vn, V^n) + {v ■ VV^n, V^n) | 
= \{V^v ■ Vn,V^n) + {Vv ■ VVn, V^n) - ^{div v, jV^nl^)] 

< \V^v\L6\Vn\L3\V^n\L2 + | VvU^o | V^nj^a 

< C(|V^t'|L2|Vn|/fi|V^n|i2 + iV^wlj/ilV^nl^^) 

< C5{\V^v\l, + \V^n\l,). 



Similarly, we have 

Putting these estimates into (2.20), we get 

^||W,^,£;)|i. + C|V^^|i.<C5(|V^^|l,, + |VV,^)li^)- (2.21) 

Combining (2.19) and (2.21) yields (2.17) if 6 is small enough. This completes the 
proof of the lemma. 

In the following lemma we give the dissipation on |Vn|iji. 

Lemma 2.3. Under the assumption (2.4), we have 

-^{div V, n) + C\Vn\l2 < C| V^;|l^i + C5\VE\l2, (2.22) 

^{div V, An) + C\VMh < C| V^;|^2 + CS\W^E\l2. (2.23) 
dt 

Proof. Notice that the condition div{pF'^) — for alH > gives 

divdiv[{l + n){E + l f] = 0, V t > 0. 

Thus we have 

=divdiv{E'^) 

= divdiv[(l +n){E + If] - divdiv{nl + E^) (2.24) 
= —An — divdiv(nE). 

Thus by applying Vj to (2.2b) and summing over i, we have 

{div v)t - (2/x + X)Adiv v + 2An = div gi, (2.25) 

where 

gi^g- div{nE). 

Multiplying the above equation by n, and then integration over IR^, we have 
I Vn|^2 = {{div v)t, n) — ((2yLt + X)Adiv v, n) — {div gi, n) 



— ^{div V, n) — {div v, tit) + ((2// + X)Av, Vn) + (^fi, Vn) 



— ■;]]:{div v, n) — {div v, f) + {div v, v ■ Vn) + \div v\2^2 
+ {i2fi + X)Av,Vn) + {gi,Vn), 

where we use the the continuity equation (2.2a). By Sobolev's, Holder's and Cauchy's 
inequalities, we obtain 

--j[{div v,n) + \Vn\l2 

< C(|Vf 1^21/1^2 + \Vv\L2\v\Ls\Vn\L3 + \Vv\l2 + \V^v\L2\Vn\L2 
+ |fi'i|L2|Vn|L2) 

< C\Wv\jj, + \\Vn\l2 + CS\W{n,E)\l2, 



which gives (2.22) if 5 is sufficiently small. 

Multiplying (2.25) by An, and then integrating over R^, we have 

2| An|^2 = -{{div v)t, An) + ((2^ + X)Adiv v, An) + {div gi, An) 
= -ftidiv V, An) + {div v, Ant) + ((S/x + \)Adiv v, An) 
+ {div gi, An) 

— —-^{div V, An) + {div v, A/) — {div v, A{v ■ Vn)) 

— {div V, Adiv v) + ((2// + X)Adiv v, An) + {div gi, An) 

— —-^{div V, An) — (SI div v, V/) — {Adiv v, v ■ Vn) 
+{Vdiv V, Vdiv v) + ((2/x + \)Adiv v, An) + {div gi, An). 

By Sobolev's, Holder's and Cauchy's inequalities, we have 

j^{div v,An) + 2\V^n\l2 

< C{\Vh,\L2\Vf\L2 + iVMLML^'l^nlLS + \V^v\l, 
+ \V'^v\L2\An\L2 + |V5fi|L2|An|L2) 

< C\Vv\l, + ^\An\l, + C6\V\n, E)\l„ 

which gives (2.23) if 5 is small enough. This completes the proof of lemma. 

In the following lemma we give the dissipation on |V(£''^ — E)\hi. 
Lemma 2.4. Under the assumption (2.4), we have 

-|(>V, E^-E) + C\W{E^ - E)\l, < C\Wv\l, + C5|V(n, E)\l,, (2.26) 

^(W, A{E^ - E)) + C\V\E^ - E)\l, < C\Vv\]j, + C5\V\n, E)|i„ (2.27) 
where W — Vu — (Vm)^ = curl u. 

Proof. Taking (2.2c)^ - (2.2c), we have 

{E'^ - E)t + W ^h^ -h-vV{E'^ - E). (2.28) 
Note the condition F^^ViF"-^ = F^^ViF"-^ for all t > 0, which means that 

VkE'^ + E^^ViE^^ = VjE'^ + E^^ViE'^, V t > 0. (2.29) 

Thus we have 

VjVkE'^ - ViVkE^^ 
= VkV.E^'' -VkViE^' 

= VkVkE'^ - VfcVfc^^* + VkiE^'^ViE'^ - E^^ViE'^) (2.30) 

-VkiE^'^ViE^' - E^'WiE^^) 
= A{E'^ - E^') + Vk{E^^ViE^^ - E^^ViE^) - Vk{E^^ViE^' - E^^WiE^^). 
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Thus by applying curl to (2.2b), we have 

Wt-i2AW + A{E'^-E) = curlg + S, (2.31) 

where the antisymmetric matrix S is defined as 

S'^ = Vk{E^^^iE'^ - E^^ViE'^) - VkiE^^^iE^' - E^'ViE^^). 

Notice that the system (2.28)-(2.31) takes a similar form as the system (2.2a)-(2.25). 
Thus after a similar argument as Lemma 2.3, (2.26) and (2.27) follows. The proof of 
lemma is completed. 

Finally, in the following lemma we give the dissipation on |V-E|j^i. 

Lemma 2.5. Under assumption (2.4). wc have 

\VE\l2<C\V{n,E^ - E)\l, (2.32) 
\V^E\l, < q V'(n, E'^ - E)\l. (2.33) 

Proof. Combining (2.24) and (2.30), we have 

/\div E = Vdivdiv E — curlcurldiv E 



= -AVn - Vdwdiv{nE) + Acurl{E - E^) + curlS. 

Thus using the property of Riesz potential, (2.4) and (2.5), we arrive at 

\divE\l, < C{\Vn\l, + |V(^^ - E)\l, + \V{nE)\l, + \EVE\l,) 
< C\V{n, E^ - E)\l, + Cd\V^E\l„ 

and 

\\/divE\l, < C{\V^n\l, + \V^iE^ - E)\l, + \V\nE)\l, + \ViEVE)\l,) 
< C\V\n, E^ - E)\l, + CS\VE\l,. 

Under the above estimate, we may deduce from (2.29) that 

\WE\l, < \divE\l, + \curlE\l, 

< C\V{n, E^ - E)\l, + Cd\VE\l, + \EVE\l, 

< C|V(n, E^ - E)\l, + C6\VE\l„ 

and 

\V^E\l2 < \VdivE\l, + \VcurlE\l, 

< C\V\n, E^ - E)\l, + CS\V^E\l, + \V{EVE)\l, 

< C|V2(n, E^ - E)\l, + CSlWEll,. 

This proves (2.32) and (2.33), and the proof lemma is completed. 



(2.34) 
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Now we are in a position to verify (2.4). Since 5 > is sufficiently small, from 
Lemma 2.1-Lemma 2.5, we can choose a constant Di > suitably large such that 



i{Di\{n, V, E)\l, + {div V, An - n) + (W, A{E^ - E) - {E'^ - E))} 
+C{\V{n,E)\]j, + \Vv\]j,)<Q, 

for any t > 0, which implies 

\{n,v,E)\l, + I {\V{n,E)\l, + \Vv\l,)<C\{no,vo,Eo)\l2, (2.35) 

since 

Di \{n,v,E)\]j^ + {div V, An-n) + {W, A{E^ - E) - {E^ - E)) ^ \{n,v, E)\jj2. 
Then (2.35) gives (2.4). Thus we prove the global existence result of Theorem 1.1. 

3 Convergence rate of the solution 

In this section wc shall prove the decay rates of the solution to finish the proof of 
Theorem 1.1. In Section 3.1, wc list some elementary conclusion on the decay-in-time 
estimates for the linearized system and a useful inequality. In Section 3.2, we shall 
first obtain the energy inequality for the derivatives of the orders from the first to the 
third, and then we show a decay-in-time estimate for the first order derivatives, where 
the error is related to the derivatives of the higher order. Finally, by combining these 
estimates we get the optimal decay rates. 

3.1. Spectral analysis and linear estimates 

We first note that the linearized system (2.2a)-(2.25) depends only on {n.div v) 
while the linearized system (2.28)-(2.31) also depends only on (W, -B^ — E). Denote 
by A* the pseudo differential operator defined by 

and let 

m = A~^div V 
be the "compressible part" of the velocity, and 

CO = A"^>V = k~'^curl V 

be the "incompressible part" of the velocity. We finally obtain 

nt + Ad^f-v-Vn, (3.1a) 
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dt - (2// + \)Ad- 2An = A'^div gi, (3.16) 

(E'^ - E)t + Auj^h'^ -h-v V{E'^ - E), (3.2a) 

ojt - nAuj - A{E'^ -E)= A-^curl g + A'^S. (3.26) 
Indeed, as the definition of d and u, and tlie relation 

V = -A~ + A'^curl uj 

involve pseudo-differential operators of degree zero, the estimates in space iJ'(M^) for 
the original function v will be the same as for (d, a;). 

Here, we just discuss the system (3.1) for example, since the system (3.2) is the 
same as system (3.1). To use the W — L'^ estimates of the linear problem for the 
nonlinear system (3.1) and system (3.2), we rewrite the solution of (3,1) as 

U{t) = K{t)Uo + f Kit- T)G{r)dT t > 0, 
Jo 

where we use the notations 

U = [n, df, Uo = [no, dof , G=[f-vVn, A'^div g^f , 

and Kit) is the solution semigroup defined by K{t) — e*-^, t > 0, with B being a 
matrix- valued differential operator given by 



B 



-A 
2A (2// + A) A 



Now we aim to analyze the differential operator B in terms of its Fourier expression 
A and to show the long time properties of the semigroup K{t). For this purpose, we 
need to consider the following linearized system 

Ut = BU. (3.3) 

Applying the Fourier transform to system (3.3), we have 

dtU=A{OU, U{0) = Uo, 
where U{t) = U{C,t) = ^U{U),^= (6,6,6) and A{0 is defined as 

"^^^^ ^ = { 2^1 -(2^ + l)|eP 

The characteristic polynomial of ^(6 is + (2// + X)k + 2|6^, which implies the 
eigenvalues are 
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The semigroup e is expressed as 



2|g|(e'^+«-e'^-*) «+e"-*-<t-e''+* _ (2j^+A)|gp(e''+*-e''-*) 

Thus the semigroup -ft'(t) has the following properties on the decay in time, which can 
be found in [10,11,14]. 

Lemma 3.1. Let A; > be an integer and 1 < I < 2. Then for any t > 0, the 
solution U{t) — {n{t),d{t)) of system (3.6) satisfies 

\V'K{t)Uo\L^ < C{l + t)-''^'''-''\\Uo\i + \V'Uo\l^) < C(l + t)-'^(^'2;^)|(n,^;)U.nH^, 

where the decay rate is measured by 

<7(/,2;A:) = ^(l-l) + |. (3.4) 

Moreover, if |no| > cq, do = for < |^| <^ 1, then there exists a positive constant C2 
such that 

Ht)\L2 > C2(l + t)-t, 
\dit)\L2 > C2(l + t)-i 

Finally, if |(no, do)\ < ICT ^or < |^| ^ 1, then there exists a positive constant C such 
that 

\in,d)it)\L2<Cil + t)-^-^{no,do)\L'2. 



We finish this subsection by listing an elementary but useful inequality [4]: 
Lemma 3.2. If ri > 1 and r2 G [0,ri], then it holds that 

/V + t - r)-'-Hl + r)-'-^ < C{n, r2)(l + t)' 

3.2. Convergence rates 

Now we will show the energy inequality as follows: 



-r2 



Lemma 3.3. Under the assumption (2.4), let {n,v,E) be the solution to the ini- 
tial value problem (2.2), then there are two positive constants C and D2 such that if 
5 > in (2.4) is small enough, it holds 

^M{t)+D2M(t) < C\V(n,v,E){t)\l2, (3.5) 
at 



13 



where the energy function M{t) defined by (3.7) is equivalent to |V(n, v, that 
is, there exists a positive constant Ci > such that 

l-\V{n,v,E){t)\l^ < M{t) < C^\V{n,v,E){t)\l^. 

Proof. Since 5 > is sufficiently small, from Lemma 2.2-Lemma 2.5, we can choose a 
constant D2 > suitably large such that 



|{D2|V(n, V, E)\l, + {div V, An) + (W, A{E' - E))} 
+Ci\V'in,E)\l, + \V^v\l,) < C6\Vin,v,E)\L2. 



Define the energy functional 

M{t)^D2\V{n,v,E%^ + {divv,An) + (W,A{E^ -E)), (3.7) 

for any t > 0, where it is noticed that M{t) is equivalent to |V(ri, v, -E)]^! since D2 can 
be large enough. Adding \V{n,v, E)\ to both sides of (3.6) gives (3.5). This completes 
the proof of the lemma. 

If we define 

N(t)= sup{1 + t)^M{t), (3.8) 

0<T<t 

then 



\Vin,v,E){t)\Hi < C^/mtj <C{l+t)'i^/Nit}. (3.9) 

To close the estimate (3.5), we shall estimate the decay rate of the first order deriva- 
tives, this will be based on Lemma 3.1 about the decay estimates on the semigroup 
K{t). Precisely, we have the following lemma. 

Lemma 3.4. Under the assumption (2.4), let {n,v,E) be the solution to the ini- 
tial value problem (2.2). Then we have 



\V{n,v,E){t)\L^<C{l + t)-^Ko + S^/N{t)), (3.10) 
where Kq = \{no,vo, Eo)\LinH^- 
Proof. Prom the Duhamel's principle, it holds that 

Thus from Lemma 3.1, we have 

\in,d)\<CKo{l + t)-'^^'''''^ + C [\l + t-T)-'^^''^''\\G{T)\L^ + \G{r)\m)dr, (3.11) 

Jo 
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where (t(1,2;1) = | by (3.4). By Holder's inequality and Sobolev's inequality, the 
nonlinear source terms can be estimated as follows: 

\G{t)\Lo. <C\{f-v Vn, (?i)Ui < C(5(|V(n, E)\l2 + \Vv\h^), 

\G{t)\Hi <C\{f-v Wn,g,)\H^ < CS\W{n,v, E)\h^ + C\Wn\HA^h\L^). 

Putting these estimates into (3.11), by (2.35), (3.9), Lemma 3.2 and Holder's inequahty, 
we arrive at 



I V(n, d){t) \l2 < CKo{l + tyl + C5 £(1 + t - r)-t (1 + r)"! y/N{^dT 
+Cj^{l + t-T)-l{l + r)-'^y^N{^\VMr)\L^dT 
< C{l + t)-l{Ko + 6,/Nit)) 



+Cv^/o(l + t - r)-t(l + T)-^VMr)\L^dr 

< Cil + t)-i{Ko + S,/N{t)) 

+C^/W)iIoa + t - r)-t(l + T)~^dT)Hj^ \VMr)\hdT)'2 

< Cil + t)-^Ko + 6./N{t)). 

Similarly, we have 

\V{uj,E^ -E){t)\L^ < C{l + t)-^Ko + S^/N{Fj). 

Combining the above two inequalities, Lemma 2.5, and the relation of v and {d,uj), we 
get (3.10). This completes the proof of the lemma. 

Now we are in a position to prove (1.6) -(1.7) in Theorem 1.1. Applying the Gron- 
wall's inequality to the Lyapunov-type inequality (3.5), by (3.10), we get 

M{t) < M{0)e-^^' + C /o e-^2(*-")|V(n, 't;,E)||,dr 

< M(0)e-°2* + C jr(Je-^^(*-^)(l + ryHxi + 5^N{r))dr 

< C(l + t)-i(M(0) + Ki + 5''N{t)). 

In view of (3.8), we have 

N{t) < C(M(0) + + S''N{t)), 

which implies 

Nit) < C(M(0) + K^) < CKl 
since > is sufficiently small. Thus (3.9) gives 

\V{n,v,E){t)\m<CKa{l + t)--^. (3.12) 

This proves (1.7). Now for (1.6), first by Sovolev's inequality and (3.12), we have 

\{n,v,E){t)\L^<C\V{n,v,E){t)\L2<CKo{l + tY-^. (3.13) 
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Meanwhile, using Lemma 3.1 and Lemma 3.2, it follows from the Duhamel's principle 
that 

|(n,d)(t)U. 

< CKoil + + C/o (1 + t - r)-||(/ - V ■ Vn,g,){T)\LinL^dT 

< CKo{l + t)~'l +C6 hi + t-T)-i\V{n,v,E){T)\Hidr (3.14) 

•J J. o c 



< CKo{l+t)-^ +C6KoJ^{l + t-T)-'^l + T)-Ur 

< CKo{l + t)-t 



Similarly, we have 

\{uj,E^-E){t)\L2<CKo(l + t)--'. (3.15) 
Finally, we derive the time-decay-rate on \E{t)\L2. From (2.29) and (2.34), we have 

\A~^curl E{t)\L2 < C\A-\EVE){t)\L2 < C\E{t)VE{t)\e < CKo{l+t)~^, (3.16) 

\A-^div E{t)\L2 < C{\{n,E-E^){t)\L2 + Ht)E{t)\L2 + \A-^S{t)\L^) 

< CKo{l + t)-'i + C\A-\EVE){t)\L2 

< CKo{l + t)-i. 

The above two inequalities give 

\E{t)\L2<CKo(l + tyi. (3.17) 
Combining (3.14), (3.15), (3.17) and the relation of v and {d,u), we have 

\{n,v,E){t)\L^<CKo{l + t)-l (3.18) 
Hence, by the interpolation, it follows from (3.13), (3.18) that for any 2 < p < 6 

|(n,^,£;)(t)Up < |(n,^,£;)(i)|l.|(n,^,E)(i)|^/ < C(l + t)-i(^4), 
where 9 — The proof of (1.6)-(1.7) is completed. 

4 Lower bound time decay rate 

In this section, we investigate the lower bound time decay for global solutions. 
Define 

g{t, x) = pit, x) — 1, m{t, x) = pu, T = pF — I. 
Then the condition divJ^ = ensures that 

div{pFF^) = diviiJ' + /)(^)] 

= divJ' -Vg + div{ -'^+^^^-f^^+<''' ), 

and 

divdivF — divdivT^ — 0. 
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Thus we have the following system which only depends on {g, div m) 



Qt + div m = 0, 



{div m)t 



{2ii + X)/\{div m) + (1 + a)/\Q = G 



1) 



where 



Gi^ adivdivi^^^^ 




By Holder's inequality and Sobolev's inequality, it is easy to verify that 



Thus by Duhamel's principle, Lemma 3.1 and the condition (1.8), we have 



\(Q,A-^div m)(t)\L2 

> \K{t){go,A-'divmo)\L^- \K{t - t){0, A-^Gi{T))\L^dT 

> c,{l + t)-i-Cj^{l+t-T)-H\A-'G,{T)\Loo + \A-'G,{T)\L2)dT 

> Ci(l + t)-t - (1 + t - T)-t (1 + T)-ldT 

> C2{l+t)--4. 



Hence (1.9) is proved. 

On the other hand, the condition F^^ViF'^^ — F'-^ViF'^'^ means that 



(|A-2Gi(t)Uoo + |A-iGi(t)|i2) < C\{g,m,J')\l, < €{1 + 1)-^ 




1+e 



where 
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Using the fact div = 0, we have 



v.(-F'^vK5^) + v,(^)) 
-v.(J-''=vKg^) + v.(^)) 

-v.(vKJ-'^g^) + v.(^)) 

+V.(V,g^)-V,(V,fg). 
Thus by applying curl to (1.1b) we have 

[curl m)t — fJ,A{curl m) + aA{J^^ — F) = Hi 

where 

= V,(^Ag) - V,(^Ag) + V,{div{pu ® u)) - V,{div{pu ® u)) 
+V.(VK-F^'=g) + V.(^)) - V.(V,(^'^S) + V,(^)) 
-V.(VK-F''=g^) + V.(^)) + + V.(-g^)) 

+V.(V.g^)-V,(V.g^). 

We also note that 

V X (u X pF'^) = pF^ ■ Vm - pF^dii' M - M ■ V(pF'^) + udw{pF'^) 
= pF^{Vuf - pF^div u- (u- V p)F^ - p{u ■ VF'^) 
= p{F^{Vuf - u ■ VF^) - {pdiv u + u- V p)F'^ 
= pFf + p,F^ = (pF^)„ 

where we used the condition div{pF'^) = 0. Thus we have 

(J^^ — J^)t + curl m — i?2, 

where 

if2 = V X (m X T'^) + cwrZ (^m) - (V x (m x T'^) + cwrZ (^m))^. 
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By Holder's inequality and Sobolev's inequality, it is easy to verify that 
(|(A^i,A^2)W|loo + |(A-^i/i,i^2)W|L0 < C-K^.m, < + 

Duhamel's principle, Lemma 3.1 and the condition (1.8), we have the following esti- 
mates of system (4.2)-(4.3): 

\{J^ - J^,K-^curi m){t)\L2 

< \K{t)iT^ - T,,A-'curl mo)\L^ - jl\K{t - T){H2{T),K-^H^{T)\L^dr 

< C + CXil+J-r)-l 

x{\{A-m,,A-^H2){T)\Lo^ + \{A-'H,,H2){T)\L2)dr ^ " ^ 

< C(l + + C/o(l + t - r)-t(l + T)--2dT 

< c(i + t)-°^^"(i+t'!). 

Combining (4.1) and (4.4) gives 

\'m-{t)\L'2 > \A^^div m(t)|x,2 — \A~^curl m{t)\L2 

> C3(l + t)-i - C{1 + t)— 

> C4(l+/)-3. 

Hence (1.10) is proved. 
By (4.4), we also have 

\{E^ - E){t)\L2 < C|(^^ - ^)(t)U. < C(l + t)— (4 5) 

Thus from (2.34), (4.1) and (4.5) , we obtain 

\A-'div E{t)\L2 > \n{t)\L2 - \n{t)E{t)\L2 - |(^^ - E){t)\L2 - \A-'S{t)\L2 

> C2(l + t)-i - C(l + t)---(i+!.!) 

> C5(l+t)-t. 

Combining the above inequality with (3.16), we arrive at 

\E{t)\L2 > \A-'^div E{t)\L2 - \A-^curl E{t)\L2 > C6(l + t)"i 
Thus, (1.11) is proved and this completes the proof of Theorem 1.1. 
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